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Abstract 

Recently, a tentative 130 GeV y-ray line signal was identified by quite a few 
groups. If correct it would constitute a "smoking gun" for dark matter annihila- 
tions. Interestingly, the spectra of the cosmic ray electrons detected by PAMELA 
and Fermi-LAT both show tiny wiggle-like structure at ~ 100 GeV, which might 
indicate a weak signal of the annihilation of ~ 130 GeV dark matter particles 
into electrons/positrons with a velocity-weighted cross section (crv) xx ^ e — 4 x 
10" 26 cm 3 s" 1 . The prospect of identifying such a potential weak dark-matter- 
annihilation electron and/or positron component by AMS-02, a mission to mea- 
sure the high energy cosmic ray spectra with unprecedented accuracy, is investi- 
gated. 

1. Introduction 

Gamma-ray line is generally thought to be a smoking gun observation of dark 
matter (DM). Recently, Bringmann et al. [Q]] and Weniger [Q] reported that there 
might be hint of a monochromatic y-ray line with energy ~ 130 GeV in the data 
recorded by Fermi Large Area Telescope (Fermi-LAT) Q3J. This y-ray line could 
be explained by ~ 130 GeV DM particle annihilation, with the velocity-weighted 
cross section (crv) xx ^ r7 ~ 10~ 27 cm 3 s _1 . This phenomenon was confirmed by a 
series of independent analyses It was argued that such a line-like struc- 

ture might originate from astrophysical emission related with the Fermi bubbles 
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|0] but the morphology analysis indicated that it is not the case 0]. Based on 
the identified spectral and spatial variations of rich structures of the diffuse y-ray 
emission in the inner Galaxy, Boyarsky et. al. argued against the DM origin of 
these structures However, the DM origin of the y-ray line emission has been 
strengthened by Su & Finkbeiner [6] and by Yang et al. The independent anal- 
yses to search for y-ray lines in the Milky Way halo by Fermi-LAT collaboration 
@] and in dwarf galaxies by [TToh found no significant signal, but the constraints 
are not tight enough to exclude such a y-ray line signal. It was also proposed that 
such a line-like signal could be tested with high energy resolution detectors in the 
near futur e 111 ill . Several models had been proposed to explain this tentative line 
structure [03^ 

Several years ago ATIC experiment discovered significant excess in the e + + e~ 
energy spectrum between 300 - 800 GeV, moreover the e + + e~ energy spectrum 



also showed possible wiggle-like structure at ~ 100 GeV 111311 . which has been 



studied by a few research groups ||14|| . The e spectra measured by PAMELA 



and Fermi-LAT both revealed tentative fine structure above ~ 100 GeV It 1 5L 1 1 6TI . 
Therefore a natural question one would ask is whether there is any connection 
between the 130 GeV line-like structure of y-ray s and the wiggle structure of 
electrons. 

In this Letter, we focus on the PAMELA/Fermi-LAT electron data and show 
that the DM scenario with mass ~ 130 GeV corresponding to the possible y-ray 
line, might contribute to the tentative fine structure of the e~ spectra around 100 



GeV. Together with the PAMELA/Fermi-LAT positron fraction [|16l Il7ll data, we 



set a constraint on the velocity-averaged cross section (crv) xx ^ e + e - < \0~ 25 cm 3 
s~ ! , consistent with all of the current bounds of the indirect detection measure- 
ments. Considering the large uncertainties of the current data, more advanced and 
dedicated experimental observations, in particular by AMS-020, are highly neces- 
sary to pin down the shape of the electron and positron spectrum, and then confirm 
or rule out a DM component in accordance with the ~ 130 GeV y-ray line. 

2. Cosmic ray propagation 



The cosmic ray (CR) propagation equation is written as follows 111811 : 
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where ip = i/r(r, p, t) is the density per unit of total particle momentum, q(r, p) 
is the source distribution function, D xx is the spatial diffusion coefficient, V = 
dV/dz x z is the convection velocity, D pp is the diffusion coefficient in momentum 
space, p = dp/dt is the momentum loss rate, T/ and j r are the time scales of 
fragmentation and radioactive decay. 

In general it is difficult to solve the propagation equation with analytical method, 
given the complicated distributions of the source, interstellar matter, radiation field 
and magnetic field. Numerical methods are developed to solve the propagation 
equations, such as GALPROP Q and DRAGON 0]. In this work we adopt 
the GALPROP package to calculate the propagation of the CR particles, includ- 
ing the contribution from DM annihilation. The diffusion-reacceleration (DR) and 
diffusion-convection (DC) models of CR propagation are adopted as illustration. 
The main parameters of propagation and source injection are compiled in Table. 
\T\ These set of propagation parameters can fit the observational B/C, 10 Be/ 9 Be 
and proton data I120ll . 



Table 1: The propagation parameters. 
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3. Model and Results 

3.1. The contribution of Dark matter annihilation and Pulsars 

From the PAMELA and Fermi-LAT data of the electrons II151 llol . we can 
see that there may be a tiny excess above ~ 100 GeV. The ATIC and Fermi- 
LAT spectra of the total electrons and positrons suggest that there is a significant 
excess at energies above 300 GeV. Therefore we adopt a three- component electron 
model to fit the data. The background electrons from primary cosmic ray sources 
contribute to the electrons below ~ 50 GeV, the pulsar component to reproduce the 
high energy excesses of the e + e~ spectra since pulsars are a kind of feasible high 
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Figure 1 : The e~ flux (left) and positron fraction (right) for DR propagation model. The dash- 
dotted (red) line is the CR background component, the dashed (green) line represents the sum of 
background and pulsar components, and the solid (blue) line is the sum of the above two and an 
additional DM component with mass ~ 130 GeV. The annihilation cross section of the DM is 
3.9 x 10~ 26 cm 3 s _1 . As shown in the right panel, due to lower flux of the corresponding "back- 
ground", it will be easier to identify a DM- orig in positron component than the electron component. 
References of the data are: PAMELA filo, Fermi-LAT IM. 



energy positron and electron source that has been widely discussed in literature 
(e-g-, PILI22ID . and the annihilation of DM particles with mass ~ 130 GeV. 
The source function of electrons and positrons from DM annihilation is 



q(E, r) 



((TV) 



xx- 
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dE 



xp 2 (r), 



(2) 



where m x is the particle mass of DM, p(r) is the spatial distribution of energy 
density, and dN/dE is the electron and positron yield spectrum produced by one 
pair of DM annihilation. In this work we use the Einasto DM density profile (12311 



p(r) = p* exp 

a* 



— I -1 



(3) 



where a* = 0.17, r* « 15.7 kpc and p* « 0.14 GeV cm -3 . 

We assume that the high energy electrons/positrons are generated through the 
cascade of electrons accelerated in the magnetosphere of pulsars 112 ll 12411 . The en- 
ergy spectrum of e + e~ injected to the galaxy from pulsars can be parameterized as 
a broken power-law with the cutoff at E c , dN/dE oc A psi E~ a exp(-E/E c ), where 
E c ranges from several tens GeV to higher than TeV, according to the models and 
parameters of the pulsars [|2l[ I25Q . And the power-law index a ranges from 1 to 
2.2 depending on the gamma-ray and radio observations D24J. In this Letter, we 
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Figure 2: Same as Figure[T]but for the DC propagation model. The DM annihilation cross section 
is 4.0 x 10~ 26 cm 3 s" 1 . 



adopt E c = 860 GeV and a = 1.28 following |27[ . 
pulsars can be parameterized the following form [|26|l 



f(R,z)cc\?L\ exp 



HP-Re) 



The spatial distribution of 
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where R Q = 8.5 kpc is the distance of solar system from the Galactic center, 
z s « 0.2 kpc is the scale height of the pulsar distribution, a = 2.35 and b = 5.56. 
The other parameter appear in the above equation is the normalization factor A psr , 
and it will be determined in our modeling. 



3.2. The Results 

Figure \T\ shows the results of the three-component model in the DR propaga- 
tion model. The CR electron background is calculated using the parameters in 
Table [Q Note for the background positron (and the secondary electron) flux we 
multiply a constant factor c e + = 1.4 to better fit the data P7Q, which may account 
for the uncertainties of the propagation model, interstellar gas distribution and the 
inelastic hadronic interaction model. To judge the improvement of the fit through 
adding the ~ 130 GeV DM component, we calculate the^ 2 value of the model. To 
minimize the effect of solar modulation, we use the data with energies higher than 
10 GeV for^- 2 calculation, although we employ a force field approximation [28] 
to approach the solar modulation with modulation potential ~ 400 MV. The^ 2 for 
the background + pulsar model (null hypothesis) is 20.9. With the presence of a 
~ 130 GeV DM component, the minimum^ 2 found is 18.9, and the best-fit cross 
section is (crv) xx ^ e + e - = 3.9 x 10~ 26 cm 3 s _1 . The fit is slightly improved but not 
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Figure 3: The 2cr upper limits on (trv)^-» e + e - for DR (dash-dotted) and DC (solid) propagation 
models respectively. 



significantly enough in case of an additional degree of freedom. The small differ- 
ence between the^ 2 inferred in these two scenarios suggests that the background 
+ pulsar model is enough to describe the data. 

Figure |2]is the same as Figure \T\ except that it is for the diffusion-convection 
propagation model. Similarly we calculate the likelihood of the DM component. 
The x 2 f° r m e null hypothesis is 21.5, while it is 19.9 in the presence of a 130 
GeV DM component. The best fit cross section is 4.0 x 10~ 26 cm 3 s _1 . As in the 
DR model, the background plus pulsar model gives reasonable fit to the data. 

It is shown above that adding a DM component do not significantly improve 
the fitting, so we turn to set an upper limit of the DM component instead. The 
95% confidence level upper limits of the annihilation cross section (crv) xx ^ e + e - 
for different mass m x are shown in Figure |3] For m x = 130 GeV, the 2<x upper 
limit of the cross section is about 8 x 10 -26 cm 3 s _1 . This constraint is stronger than 
that derived through y-ray observations of the Milky Way halo One caution 
is that the constraints also depend on the assumptions used for the background, 
which are not same as in [Ml. 



Since the current observational data are not precise enough to see whether 
there is a tiny structure of the electron spectra at ~ 100 GeV, we would like to 
discuss the potential of the AMS-02 experiment on this issue. We use Monte 
Carlo simulation to produce the expected electron spectra of AMS-02, based on 
the theoretical fluxes of electrons (f> e for the background + pulsar + DM scenario 
as given in Figures [T] and |U The events number of AMS-02 can be estimated as 
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L29] 



N = AtA e f dE fdE' <p e (E') — ^e"^, (5) 
Jae J \2jrcr 

where cr = ^(0.106/ V^/GeV) 2 + 0.0125 2 x£/2B represent the energy res- 
olution, At is the operating time, (p e is the electron flux, and A e is the geometrical 
acceptance of electron which is taken to be 0.045 m 2 sr 113 IN . The electrons are 
binned into 25 bins logarithmically from 10 GeV to 1 TeV. The observed number 
of events in each energy bin is generated based on a Possion distribution with ex- 
pected value N calculated by Eq. (5). The logarithm of the likelihood function is 
defined as 

ln£ = ^AUn^-^-lniV,-!, (6) 

where Nj is the simulated observational electron counts in energy bin i, and <p\ is 
the expected counts in the energy bin. The test statistic of the DM signal is defined 
as TS = -2 ln(£nuii/-£bestX where Xnuii is the best fit likelihood of null hypothesis 
(background + pulsar scenario), and Xbest is the best fit likelihood of the model 
with 130 GeV DM component. For a series of annihilation cross sections, we cal- 
culate the TS value and find the required exposure time T which makes TS ~ 25 
(approximately corresponding to 5<x significance). The results are summarized 
in Table [2] It is shown that if the DM annihilation cross section is larger than 
10" 26 cm 3 s -1 , this 130 GeV DM signal may be identified by AMS-02. 



Table 2: The predicted operating time for AMS-02 to identify the 130 GeV DM annihilating into 
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We may further note that the DM signal, if exists, will be more prominent in 
the positron fraction than in the electron spectrum. As shown in Figures [TJ and 
[H the relative excess of the DM component is only ~ 2% compared with the 
background + pulsar flux. However, for the positron fraction the contribution of 
the DM component could be more than 10% (relative to the background + pulsar 
components) for the best fit cross section 4 x 10 26 cm 3 s _1 . 
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3.3. Constraints from other observations 

The observations of y-rays (the internal bremsstrahlung and inverse Compton 
radiation component) and/or radio emission (the synchrotron radiation compo- 
nent) may constrain the current scenario that DM annihilates into electrons/positrons 
n32l l9fl. The lack of evident continual spectrum component associated with the 130 
GeV y-ray line from the Galactic center suggests that the DM particles mainly 
annihilate into final states with few y-rays (such as e + e~, n + u~, or neutrinos), oth- 
erwise the observed relic density can not be explained (13311 . Such a speculation 
is one of the main motivation of this work. It was shown that the latest y-ray ob- 
servations by Fermi-LAT can constrain the (<xv)^_> e + e - to the level of 10~ 24 cm 3 
s for DM mass ~ 100 GeV §], which is much larger than the constraint (i.e., 



(crvL^ e + e - < 10 25 cm 3 s l ) yielded in this work. 



4. Conclusion 

The spectra of the CR electrons detected by PAMELA and Fermi-LAT both 
show small wiggle-like structure at ~ 100 GeV, potentially consisting of a weak 
signal of ~ 130 GeV DM particles annihilating into electrons/positrons with a 
cross section (crv) xx ^ e + e - < 10" 25 cm 3 s _1 . It maybe connect with the recently 
reported ~ 130 GeV y-ray line emission which may be the result of DM annihi- 
lation in the Milky Way. We investigate the contribution to the electron spectrum 
and positron fraction of such a DM component. It is found that adding the 130 
GeV DM component can improve the fit to the data, but the improvement is not 
significant enough. We further use the current data to set constraints on the DM 
annihilation cross section to e + e~. 

As found in our modeling, the DM-origin electrons/positrons, if there are, 
are likely less than ~ 2% (4%) of the background + pulsar e~ (e~ + e + ) flux at 
~ 130 GeV. Hence, accurate measurements of the spectrum of the electrons (and 
positrons) by AMS-02 and DAMPE/CALEli are highly necessary to test the DM 
origin of the wiggle-like structure (see Table 2 for the expected performance of 
AMS-02 in order to identify such a component). Finally we would like to point 
out that it may be easier to identify a DM-origin positron component than the 
electron component due to lower flux of the corresponding "background". 



3 DAMPE and CALET (http://calet.phys.lsu.edu/) are able to measure the total spectrum of 
cosmic ray electrons and positrons accurately. 
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